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Abstract.  This  paper  presents  an  evaluation  and  validation  of 
the  Naval  Research  Laboratory’s  COAMPSw  real-time  fore¬ 
casts  during  the  VOCALS-REx  over  the  area  off  the  west 
coast  of  Chile/Peru  in  the  Southeast  Pacific  during  Octo¬ 
ber  and  November  2008.  The  analyses  focus  on  the  marine 
boundary  layer  (MBL)  structure.  These  forecasts  are  com¬ 
pared  with  lower  troposphere  soundings,  in  situ  surface  mea¬ 
surements,  and  satellite  observations.  The  predicted  mean 
MBL  cloud  and  surface  wind  spatial  distributions  are  in  good 
agreement  with  the  satellite  observations.  The  large-scale 
longitudinal  variation  of  the  MBL  structure  along  20°  S  is 
captured  by  the  forecasts.  That  is,  the  MBL  height  increases 
westward  toward  the  open  ocean,  the  moisture  just  above  the 
inversion  decreases,  and  the  MBL  structure  becomes  more 
decoupled  offshore.  The  observed  strong  wind  shear  across 
the  cloud-top  inversion  near  20°  S  was  correctly  predicted  by 
the  model.  The  model’s  cloud  spatial  and  temporal  distribu¬ 
tion  in  the  15  km  grid  mesh  is  sporadic  compared  to  satellite 
observations.  Our  results  suggest  that  this  is  caused  by  grid- 
scale  convection  likely  due  to  a  lack  of  a  shallow  cumulus 
convection  parameterization  in  the  model.  Both  observations 
and  model  forecasts  show  wind  speed  maxima  near  the  top 
of  MBL  along  20°  S,  which  is  consistent  with  the  westward 
upslope  of  the  MBL  heights  based  on  the  thermal  wind  rela¬ 
tionship.  The  forecasts  produced  well-defined  diurnal  varia¬ 
tions  in  the  spatially-averaged  MBL  structure,  although  the 
overall  signal  is  weaker  than  those  derived  from  the  in  situ 
measurements  and  satellite  data.  The  MBL  heights  are  gen- 
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erally  underpredicted  in  the  nearshore  area.  An  analysis  of 
the  sensitivity  of  the  MBL  height  to  horizontal  and  vertical 
grid  resolution  suggests  that  the  underprediction  is  likely  as¬ 
sociated  with  overprediction  of  the  mesoscale  downward  mo¬ 
tion  and  cold  advection  near  the  coast. 


1  Introduction 

The  marine  boundary  layer  is  known  to  play  key  roles  in  reg¬ 
ulating  large-  and  meso-scale  atmospheric  circulations  as  it 

serves  as  the  major  source  of  moist  available  energy  for  at¬ 
mospheric  circulation  and  the  major  sink  of  atmospheric  ki¬ 
netic  energy  (e.g.,  Lorenz,  1978;  Palmen  and  Newton,  1969). 
Clouds  are  frequently  present  within  the  MBL  in  the  form  of 
overcast  stratocumulus/stratus,  organized  scattered  shallow 
cumulus  clouds  or  transitional  mixture  of  both;  they  cover  a 
third  of  the  world’s  oceans  and  form  a  persistent  feature  of 
the  Earth’s  climate  and  regional  weather  (Klein  and  Hart¬ 
mann,  1993;  Felsch  and  Whitlatch,  1993).  These  clouds 
have  a  significant  influence  on  the  underlying  surface  en¬ 
ergy  budget  by  significantly  reflecting  incoming  solar  energy, 
but  only  moderately  increasing  downward  longwave  radia¬ 
tive  flux  (Klein  and  Hartmann,  1993).  In  addition,  the  cloud- 
top  radiative  cooling  created  by  the  clouds  and  the  latent  heat 
released  by  convective  updrafts  fundamentally  change  the 
nature  of  the  turbulence  dynamics  compared  to  a  cloud-free 
MBL,  leading  to  changes  in  thermodynamic  and  wind  pro¬ 

files  in  the  MBL  (Lilly,  1968).  Consequently,  air-sea  inter¬ 
action  processes  are  profoundly  impacted  by  the  presence  of 
these  clouds  at  both  global  and  regional  scales.  It  is  of  great 
importance,  therefore,  that  coupled  climate  and  mesoscale 
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models  are  able  to  realistically  predict  the  MBL  structure, 
cloud  coverage  and  their  variability  (Ma  et  al.,  1996;  Bony 
and  Dufresne,  2005;  de  Szoeke  et  al.,  2006,  and  Shulman  et 
al.,  2007). 

The  cloud-topped  MBL  is,  however,  notoriously  difficult 
to  simulate,  because  turbulence-cloud  interaction  processes 
cannot  be  explicitly  resolved  at  current  grid  resolutions  of 
global  and  regional  models  and  thus  must  be  parameterized. 
Lor  example,  an  interactive  turbulence- shallow  convection 
parameterization  is  essential  to  predict  the  transition  from 
stratocumulus  to  shallow  cumulus  clouds  (Park  and  Brether- 
ton,  2009).  Lurthermore,  mesoscale  cloud  variability  may 
also  be  driven  by  the  so-called  aerosol-cloud-precipitation  in¬ 
teraction  mechanism  proposed  by  Albrecht  (1989),  in  which 
low  CCN  (Cloud  Condensation  Nuclei)  concentration  pro¬ 
motes  more  precipitation,  leading  to  the  destruction  and 
structural  change  of  the  clouds  (e.g.,  Stevens  et  al.,  2005; 
Wang  and  Leingold,  2009).  The  representation  of  this  com¬ 
plex  interaction  process  has  not  yet  been  developed  for  use 
in  regional  and  global  models. 

Two  primary  focuses  of  VOCALS-REx  (Variability  of 
the  American  Monsoons  Ocean-Cloud-Atmosphere-Land 
Study-Regional  experiment)  are:  (1)  the  aerosol-cloud- 
precipitation  interaction  in  the  MBL;  and  (2)  chemical  and 
physical  couplings  between  the  upper  ocean,  the  land,  and 
the  atmosphere  (Wood  et  al.,  2010).  The  field  campaign  took 
place  off  the  west  coast  of  Chile/Peru  in  the  Southeast  Pacific 
during  October  and  November  2008.  Comprehensive  mea¬ 
surements  were  made  on  many  platforms  including  aircrafts, 
ships,  satellites,  buoys,  and  land-based  observation  sites  to 
provide  description  of  the  MBL  structure,  turbulence  dynam¬ 
ics,  aerosol-cloud  microphysics  and  air-sea  interaction.  Par¬ 
ticularly,  intensive  observations  were  carried  out  along  the 
latitude  20°  S  to  provide  detailed  time-height  and  longitude- 
height  cross-sections  of  the  MBL. 

One  important  motivation  for  VOCALS-REx  is  to  improve 
the  prediction  of  MBL  clouds  by  both  global  and  regional 
models.  Because  of  the  prominent  impacts  of  MBL  clouds  on 
global  climate  as  discussed  earlier,  modeling  efforts  have  tra¬ 
ditionally  emphasized  global  models.  In  the  last  few  decades, 
however,  an  increasing  number  of  studies  have  focused  on 
the  regional  modeling  investigations  of  MBL  clouds,  in¬ 
cluding  regional  simulations  using  a  simple  two-layer  model 
(Wang  et  al.,  1993),  cloud  evaluation  and  process  studies  us¬ 
ing  mesoscale  models  in  a  regional  climate  mode  (e.g.,  Wang 
et  al.,  2004a,  b;  McCaa  and  Bretherton,  2004;  Garreaud 
and  Munoz,  2004),  investigations  of  cloud  and  dynamic  pro¬ 
cesses  in  case  studies  (Mocko  and  Cotton,  1995;  Mechem 
and  Kogan,  2003;  Thompson  et  al.,  2005),  and  study  of  shal¬ 
low  cumulus  effects  in  regional  ocean- atmosphere  coupled 
climate  model  (de  Szoeke  et  al.,  2006).  These  regional  model 
studies  demonstrate  that  the  fidelity  of  the  MBL  cloud  sim¬ 
ulations  is  highly  dependent  on  the  quality  of  the  cloud  and 
turbulence  parameterizations  as  well  as  the  accuracy  of  the 
large-  and  mesoscale  flows.  A  recent  assessment  of  current 


global  and  regional  models  clearly  shows  weakness  in  the 
simulated  cloudiness  and  MBL  heights,  two  important  pa¬ 
rameters  in  the  MBL  cloud  prediction  (Wyant  et  al.,  2010). 

During  October-November  2008,  NRL’s  (Naval  Research 
Laboratory’s)  COAMPS  (Coupled  Ocean/Atmosphere 
Mesoscale  Prediction  System)  real-time  forecasts  were 
performed  in  support  of  VOCALS-REx.  The  current  work 
presented  here  focuses  on  evaluation  and  validation  of 
these  MBL  forecasts  using  observations  taken  during  the 
VOCALS-REx  field  program.  It  is  motivated  by  the  follow¬ 
ing  questions:  to  what  degree  do  the  predicted  cloud  fields 
such  as  integrated  cloud  water  path  agree  with  observations? 
To  what  degree  does  the  predicted  MBL  structure  agree  with 
that  observed?  What  specific  aspects  of  the  model  can  be 
improved  to  better  simulate  the  MBL  cloud  system? 

Section  2  describes  COAMPS,  the  forecast  setup,  and  the 
observational  datasets.  Sections  3-6  analyze  COAMPS  re¬ 
gional  fields,  MBL  variations  along  20°  S,  statistics  of  the 
surface  variables,  and  mean  MBL  diurnal  structure.  Sec¬ 
tion  7  discusses  the  issue  of  the  model  grid  resolution  and 
Sect.  8  gives  the  summary  and  conclusion. 

2  COAMPS  real-time  forecast  setup  and  observation 
datasets 

The  atmospheric  component  of  COAMPS  represents  both  a 
forecast  model  and  a  complete  three-dimensional  data  assim¬ 
ilation  system.  The  forecast  model  is  a  finite-difference  ap¬ 
proximation  to  the  fully  incompressible  and  non-hydro  static 
equations  that  govern  atmospheric  motions  (Hodur,  1997). 
The  main  model  physical  representations  relevant  to  this 
study  are  briefly  described  as  follows.  Lor  the  boundary 
layer  prediction,  the  surface  layer  is  represented  using  the 
modified  version  of  the  Louis  (1982)  surface  flux  parame¬ 
terization  as  formulated  by  Wang  et  al.  (2002);  turbulence 
mixing  processes  are  parameterized  following  the  Mellor- 
Yamada’s  1.5  order  turbulence  closure  model  that  predicts 
turbulence  kinetic  energy  and  determines  the  turbulent  fluxes 
in  terms  of  the  down-gradient  transport  approach  (Mellor 
and  Yamada,  1982).  To  account  for  latent  heat  effect  of 
clouds  on  the  buoyancy  in  turbulence  mixing,  the  subgrid- 
scale  cloud  fraction  is  calculated  using  a  Gaussian  distribu¬ 
tion  of  turbulent  fluctuations  of  conserved  variables  (Burk 
and  Thompson,  1989).  However,  the  cloud  fraction  is  not 
used  in  the  radiation  calculation.  The  cumulus  parameter¬ 
ization  of  Kain  and  Lritsch  (1990)  is  used  to  represent  the 
convection  at  the  subgrid  scale;  and  it  is  only  activated  when 
the  grid  size  is  greater  than  10  km.  The  current  version  of 
this  convection  scheme,  however,  is  not  designed  to  represent 
shallow  cumulus  clouds,  which  is  a  limitation  that  will  be 
discussed  further  in  Sect.  4.  Clouds  and  precipitation  at  grid 
scale  are  predicted  using  the  single  moment  scheme  of  cloud 
microphysics  formulated  by  Rutledge  and  Hobbs  (1983). 
Recently,  the  four-stream  radiation  parameterization  of  Lu 
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and  Liou  (1992)  was  implemented  in  COAMPS  (Liu  et  al., 
2009). 

The  initialization  of  COAMPS  forecasts  is  performed  by  a 
data  assimilation  system  that  uses  multivariate  optimum  in¬ 
terpolation  analyses  of  soundings,  surface  and  satellite  data 
and  blended  with  the  previous  12-h  COAMPS  forecast  fields. 
Lateral  boundary  conditions  for  the  outermost  domain  are 
provided  by  the  Navy  Operational  Global  Analysis  and  Pre¬ 
diction  System  (NOGAPS)  forecast  fields.  The  sea  surface 
temperature  (SST)  is  generated  by  the  NRL  Coupled  Ocean 
Data  Assimilation  System  (NCODA),  which  applies  a  three- 
dimensional,  multivariate,  optimum-interpolation  method 
and  integrates  all  available  ocean  observations  in  real  time, 
including  ship,  buoy,  and  satellite  observations  as  discussed 
in  Cummings  (2005). 

During  VOCALS-REx,  COAMPS  provided  48  h  twice- 
daily  forecasts  with  a  3-h  interval  output,  initialized  at 
00:00  UTC  and  12:00  UTC,  between  20  October  and  29 
November  2008.  As  shown  in  Pig.  1,  the  forecast  domain 
is  configured  to  contain  three  nested  grid  meshes  (45  km, 
15  km,  and  5  km)  with  45  levels  in  the  vertical,  of  which 
26  levels  reside  below  2.5  km.  Our  analysis  is  focused  on 
the  forecast  from  the  15  km  domain,  as  the  15  km  grid  mesh 
covers  most  of  the  VOCALS-REx  observation  areas  and  has 
a  relatively  high  horizontal  resolution.  Some  results  from  the 
45  km  and  5  km  grid  meshes  will  also  be  presented  to  provide 
comparisons  in  Sect.  2  (Pig.  2)  and  Sect.  7  (Pigs.  15  and  16). 

The  observational  data  employed  in  this  study  include  both 
satellite  observations  and  in-situ  measurements  made  dur¬ 
ing  VOCALS-REx.  Satellite  observations  of  liquid  water 
path  (LWP)  and  surface  wind  speed  are  from  4  passive  mi¬ 
crowave  radiometers  onboard  the  AMSR-E  (Advanced  Mi¬ 
crowave  Scanning  Radiometer  -  Earth  Observing  System), 
TRMM  (Tropical  Rainfall  Measuring  Mission)  Microwave 
Imager  (TMI),  and  the  SSM/I  (Special  Sensor  Microwave 
Imager)  P13  and  P15  satellites.  SST  fields  from  the  AMSR- 
E  and  TMI  satellites  were  used  as  an  independent  evalua¬ 
tion  of  the  NCODA  SST  fields  used  as  the  lower  boundary 
condition  in  these  COAMPS  forecasts.  Additionally,  vector 
surface  winds  from  QuickSCAT  (Quick  Scatterometer)  are 
also  used  for  comparison.  Atmospheric  sounding  data  used 
in  this  work  were  collected  on  board  NOAA  R/V  Ronald 
H.  Brown  (RHB)  along  20°  S  and  on  board  CIRPAS  (Cen¬ 
ter  for  Interdisciplinary  Remotely  Piloted  Aircraft  Studies) 
Twin  Otter  research  flights  at  the  location  of  20°  S  and  72°  W. 
Ship-based  surface  observations  and  rawinsondes  are  avail¬ 
able  every  10  min  from  the  multi-year  NOAA  tropical  east¬ 
ern  Pacific  synthesis  data  set  (de  Szoeke  et  al.,  2010).  These 
surface  observations  consist  of  a  10-min  synthesis  data  set 
from  the  VOCALS-REx  RHB  cruise  and  data  collected  at 
the  WHOI  buoy  located  at  85°  W  and  20°  S.  Readers  are  re¬ 
ferred  to  Wood  et  al.  (2010)  for  a  complete  description  of  the 
in  situ  observations  and  platforms  during  VOCALS-REx. 


Fig.  1.  COAMPS  forecast  domains  for  the  VOCALS-REx.  Three 
grid  meshes  (45  km,  15  km,  5  km)  have  151  x  151,  199  x  188,  and 
151  x  151  grid  points,  respectively. 


3  Regional  fields 

In  this  section,  we  first  examine  the  mean  large-scale  fea¬ 
tures  by  comparing  the  COAMPS  41 -day  (20  October  to 
29  November)  averaged  fields  on  the  45  km  grid  mesh  with 
those  derived  from  the  satellite  data  over  an  expanded  area  in 
the  Southeast  Pacific.  We  then  focus  on  the  mean  and  diur¬ 
nal  variability  of  the  wind  and  cloud  fields  using  the  models 
results  on  the  15  km  grid  mesh  over  the  VOCALS-REx  area. 

As  shown  in  Fig.  2a  and  b,  the  model  SST  accurately  cap¬ 
tures  all  the  major  features  shown  by  the  AMSR-E  and  TMI 
SST  field.  Both  SST  fields  are  colder  close  to  the  Chilean 
coast,  driven  by  the  low-level  jet  of  strong  southerly  winds 
parallel  to  the  coast  shown  in  the  wind  fields  between  40°  S 
and  25°  S  (Fig.  2c  and  d).  The  locally  warm  SST  area  cen¬ 
tered  at  72°  W  and  20°  S  results  from  the  weak  southerly 
winds  blocked  by  the  Andes.  A  broad  region  of  cold  SST 
(green  shaded)  is  present  to  the  northwest  of  the  coastal  up- 
welling  zone;  it  appears  to  result  from  the  oceanic  cold  ad- 
vection  driven  by  the  southeast  wind  from  the  low-level  jet 
area  to  the  northern  side  of  the  high  pressure  system.  Even 
though  the  model  SST  is  a  product  of  the  ocean  data  assimi¬ 
lation  system,  the  AMSR-E  and  TMI  data  were  not  included 
in  NCODA  for  the  ocean  analysis.  Therefore,  this  compari¬ 
son  also  serves  as  an  independent  validation  of  the  NCODA 
product. 

The  forecast  surface  wind  vectors  agree  with  the  satel¬ 
lite  QuickSCAT  data  in  the  following  main  features.  The 
southerly  coastal  jet  along  the  Chilean  coast  is  consistent 


www.atmos-chem-phys.net/1 1/421/201 1/ 


Atmos.  Chem.  Phys.,  11,  421-437,  2011 


424 


S.  Wang  et  al.:  Real-time  forecast  of  MBL  during  VOCALS-REx 


c)  COAMPS  Surface  Wind  Speed  d)  QuikSCAT  Surface  Wind  Speed 


e)  COAMPS  Sea-Level  Pressure  f)  COAMPS  Wind  Speed  at  1 .8-km  Height 


40°S105°W  90°W  75°w 


Fig.  2.  Mean  large-scale  fields  derived  from  the  COAMPS  forecasts 
(45  km  grid)  and  the  satellite  data,  (a)  COAMPS  SST;  (b)  AMSR-E 
and  TMI  SST;  (c)  COAMPS  10-m  wind  speed  (color  shading)  and 
vectors;  (d)  QuickScat  wind  speed  (color  shading)  and  vectors;  (e) 
COAMPS  sea-level  pressure;  and  (f)  COAMPS  1.8  km  wind  speed 
(color  shading)  and  vectors. 


with  the  satellite  winds  in  the  overall  correct  location  and 
direction  of  the  jet  axis,  although  the  simulated  jet  appears 
stronger  and  extend  further  north.  Weak  wind  speed  areas  are 
present  near  the  subtropical  high  pressure  center  at  105°  W 
and  31°  S  and  nearshore  at  20°  S;  and  the  relatively  high  wind 
speeds  exist  in  the  northern  side  of  the  high  pressure  sys¬ 
tem  (Fig.  2e).  The  Chilean  low-level  jet  often  takes  place  in 
the  austral  spring  and  summer  seasons  and  appears  to  be  as¬ 
sociated  with  the  northerly  directed  pressure  gradient  force 
along  the  coast  induced  by  mid-latitude  synoptic  conditions 
(Garreaud  and  Munoz,  2005;  Munoz  and  Garreaud,  2005). 
A  recent  study  by  Jiang  et  al.  (2010)  using  these  COAMPS 
forecasts  shows  that  the  subtropical  high  pressure  system  is  a 
major  forcing  regulating  the  Chilean  low-level  jet  variability. 

As  shown  in  Fig.  2f,  the  main  features  of  the  predicted 
winds  just  above  the  MBL  inversion  at  1.8  km  height  include 
the  anti-cyclonic  circulation  associated  with  the  subtropical 
high  pressure  system  and  the  weak  wind  speed  area  between 
the  high  pressure  system  and  the  coast.  In  particular,  a  nar¬ 
row  northerly  jet  flow  extends  from  15°  S  to  30°  S  along 
the  coast,  contrasting  the  southerly  flow  within  the  MBL 


(Fig.  2d)  and  resulting  in  a  sharp  wind  shear  across  the  in¬ 
version.  The  existence  of  the  strong  wind  shear  is  confirmed 
by  the  aircraft  observations  as  will  be  presented  later. 

We  next  examine  the  cloud  and  wind  fields  with  a  focus 
on  the  VOCALS-REx  area  using  the  15  km  grid  forecasts. 
The  satellite  data  used  for  this  comparison  are  derived  from 
the  microwave  measurements  taken  on  AMSR-E,  SSM/I  and 
TMI  satellites,  which  together  have  a  sampling  frequency  of 
4-8  times  daily  at  each  grid  point  in  the  domain  depending 
on  the  detailed  sampling  characteristics  of  each  satellite. 

Figures  3  a  and  b  show  that  the  general  patterns  of  the  mean 
liquid  water  path  (LWP)  fields  from  the  forecast  and  satellite 
data  are  very  similar.  Both  LWP  fields  are  characterized  by 
extensive  coverage  over  the  southeast  Pacific  with  relatively 
thick  LWP  (~100  to  120  g  m-2)  covering  the  western  part  of 
the  domain  and  thinner  clouds  along  the  coast.  The  LWP  de¬ 
creases  significantly  to  the  south  of  25°  S,  where  subsidence 
tends  to  dominate  near  the  center  of  the  high  pressure  system. 
The  COAMPS  underpredicts  nearshore  clouds,  particularly 
along  the  coast  between  16°  S  and  20°  S  and  to  the  south  of 
25°  S.  This  deficiency  is  likely  caused  by  MBL  heights  that 
are  too  low  compared  with  those  observed  as  discussed  in  the 
next  sections. 

As  shown  in  Fig.  3c-d,  surface  wind  speeds  from  both 
the  COAMPS  simulations  and  the  satellite  data  indicate  a 
weak  wind  area  near  the  coast  between  25°  S  and  17°  S.  Both 
tend  to  have  maximum  values  in  the  low-level  jet  region  off 
the  coast  of  Chile  and  to  the  south  of  25°  S,  although  the 
maxima  from  the  forecast  (~8ms_1)  are  stronger  than  those 
from  the  satellite  observations  (~7  ms-1).  The  center  of  the 
Southeast  Pacific  high  pressure  systems  is  characterized  by 
weak  winds  and  low  LWP  values  associated  with  strong  sub¬ 
sidence,  as  shown  in  both  the  modeled  results  and  satellite 
observations.  The  comparison  for  the  15  km  mesh  is  gener¬ 
ally  consistent  with  that  for  the  45  km  mesh. 

Diurnal  variability  in  the  cloud  field  is  a  dominant  fea¬ 
ture  over  the  southeast  Pacific  (e.g.,  Rozendaal  et  al.,  1995; 
O’Dell  et  al.,  2008;  Wood  et  al.,  2009).  It  is  primarily  driven 
by  the  cloud  solar  absorption  during  daytime,  which  tends  to 
stabilize  the  MBL  and  decouple  the  cloud  from  the  surface- 
based  subcloud  layer  (Nicholls,  1984).  To  examine  the  statis¬ 
tical  behavior  of  diurnal  variation  of  the  clouds,  a  harmonic 
regression  analysis  with  the  diurnal  cycle  is  performed  on 
both  the  forecasts  and  the  satellite  data  (Figs.  4  and  5).  Be¬ 
cause  the  sampling  rate  of  the  satellite  data  is  about  4-8  times 
per  day,  we  only  focus  on  the  dominant  24  h  diurnal  cycle. 
The  amplitude  of  the  LWP  diurnal  cycle  is  about  half  of  the 
mean  LWP  in  both  the  model  and  satellite  observations,  im¬ 
plying  the  presence  of  a  significant  cloud  diurnal  variabil¬ 
ity.  The  diurnal  amplitude  of  the  cloud  LWP  is  smaller  for 
the  forecasts  (~30-50gm-2)  than  for  the  satellite  (~40- 
60  gm-2).  The  diurnal  signals  are  noticeably  weaker  near 
the  Peruvian  coast  in  both  the  satellite  and  forecasts,  where 
the  amplitudes  from  both  fields  are  about  10-20  g  m- 2 .  The 
satellite  LWP  diurnal  amplitude  is  relatively  smooth  with 
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Fig.  3.  Comparison  of  the  mean  LWP  and  wind  speed  between 
the  COAMPS  forecasts  (15  km  grid)  and  the  satellite  data,  (a) 
COAMPS  LWP;  (b)  satellite  LWP;  (c)  COAMPS  surface  wind 
speed;  and  (d)  satellite  wind  speed. 


the  maximum  values  centered  at  81°  W  and  21°  S,  while  the 
forecast  is  rather  spotty  with  high  value  areas  being  approx¬ 
imately  parallel  to  wind  vectors.  Further  examination  of  the 
instantaneous  cloud  fields  reveals  that  this  discontinuity  is 
likely  caused  by  grid- scale  convection  and  probably  due  to 
the  lack  of  shallow  cumulus  convection  parameterization  in 
the  model,  which  will  be  discussed  further  in  Sect.  4. 

To  estimate  the  diurnal  contribution  to  the  total  LWP  vari¬ 
ability,  we  calculate  the  correlation  of  the  LWP  diurnal  har¬ 
monic  regression  with  the  satellite  and  the  forecasts,  respec¬ 
tively.  The  correlation  values  for  the  COAMPS  (~0.3-0.6) 
are  smaller  than  those  for  the  satellite  data  (^0.6-0. 8)  over 
large  areas  to  the  west  of  72°  W  and  the  north  of  25°  S.  These 
values  indicate  that  while  the  diurnal  cycle  can  only  account 
for  10-36%  of  the  total  variance  in  the  predicted  LWP,  it  ac¬ 
counts  for  as  high  as  36-64%  in  the  satellite  data  over  these 
areas.  Over  the  coastal  areas  near  20°  S,  the  modeled  corre¬ 
lation  is  about  0.6,  slightly  higher  than  the  observation.  The 
phase  of  the  diurnal  variation  in  the  forecasts  is  in  general 
agreement  with  the  satellite  data  over  large  region  to  the  east 
of  85°  W  (Fig.  4e  and  f),  although  the  timing  of  the  maxima 
consistently  appears  2-4  h  earlier  for  the  model  than  for  the 
satellite. 

Surface  wind  speeds  from  both  the  forecast  and  the  satel¬ 
lite  data  show  diurnal  signals  over  large  areas  off  the  Chilean 
coast  reaching  as  far  as  82°  W  as  shown  in  Fig.  5.  For  the 
forecast,  the  area  with  the  largest  diurnal  amplitude  (~1.2- 
1.8  ms-1)  appears  to  be  in  regions  of  relatively  weak  mean 
wind.  The  western  part  of  the  area  overlaps  the  region  of 
the  satellite’s  large  amplitudes  centered  on  (80°  W,  26°  S), 


-90  -80  -70-90  -80  -70 

Longitude  Longitude 


Fig.  4.  Comparison  of  the  harmonic  analysis  of  LWP  between  the 
COAMPS  forecasts  (left  column)  and  the  satellite  data  (right  col¬ 
umn).  (a)  Predicted  diurnal  amplitude;  (b)  satellite  diurnal  ampli¬ 
tude;  (c)  correlation  with  the  predicted  harmonic  diurnal  function; 
(d)  correlation  with  the  satellite  diurnal  function;  (e)  local  hours  of 
the  LWP  maximum  in  the  model;  and  (f)  local  hours  of  the  LWP 
maximum  in  the  satellite  data. 

whose  values  (1.2  m  s_1)  are  less  than  the  forecasts.  Both  the 
forecasts  and  the  satellite  place  the  large  diurnal  wind  speed 
amplitude  areas  to  the  north  and  northwest  of  the  low-level 
jet  core  offshore  of  central  Chile;  it  is  particularly  true  for 
model  results.  This  is  consistent  with  Garreaud  and  Munoz 
(2005)  who  conclude  that  the  diurnal  cycle  in  climatological 
region  of  the  maximum  jet  wind  speed  is  small.  The  correla¬ 
tion  pattern  of  the  wind  speed  and  its  diurnal  variation  shows 
a  general  agreement  between  the  satellite  and  the  forecasts 
with  larger  values  covering  the  coastal  area  and  extending 
offshore  (Fig.  5c  and  d).  For  both  the  COAMPS  and  satellite 
wind  speed  fields,  the  maximum  wind  speeds  occur  in  the 
afternoon  along  the  coast  between  25°  S  and  17°  S  (Fig.  5e- 
f),  which  is  likely  produced  by  the  afternoon  strengthening 
of  the  southerly  winds  (Rutllant,  1993).  For  other  areas,  the 
wind  speeds  are  maximum  during  the  late  evening  and  early 
morning.  The  diurnal  variation  of  the  surface  wind  speed 
only  accounts  for  10-25%  of  the  variances  for  the  satellite  as 
well  as  the  forecast  results,  a  percentage  much  less  than  that 
for  the  clouds. 
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Fig.  5.  The  same  as  in  Fig.  4  except  for  the  surface  wind  speed  and 
the  correlation  coefficient  ranges  0-0.5  now.  Note  that  the  scale 
in  the  color  coding  used  for  the  correlation  coefficients  is  different 
from  the  corresponding  one  used  in  Fig.  3c  and  d. 


There  are  some  apparent  differences  between  the  forecast 
and  satellite  wind  fields  along  west  coast  of  Peru.  The  45-km 
grid  nest  wind  speed  reaches  a  local  maximum  (~9ms_1) 
compared  to  only  7  ms-1  from  the  QuickSCAT  data  (Fig.  2), 
although  the  15 -km  grid  nest  forecast  compares  better  with 
the  combined  satellite  data  from  AMSR-E,  SSM/I  and  TMI 
satellites.  This  may  be  partly  due  to  the  fact  that  QuikSCAT 
only  has  1-2  passes  over  a  location  each  day,  while  the  com¬ 
bined  dataset  has  4-8  passes.  The  forecast  diurnal  amplitude 
over  the  area  off  Chile  between  76°  W  and  72°  W  is  ~1.25- 
1.75  ms-1,  significantly  larger  than  that  from  the  satellite 
data,  ~0.5ms_1.  Previous  multi-year  satellite  data  analy¬ 
ses,  however,  also  show  that  the  morning-to-afternoon  wind 
speed  increase  can  be  as  large  as  1.5 ms-1  in  the  same  area 
(Munoz,  2008).  The  diurnal  phase  difference  between  the 
forecast  and  the  satellite  is  about  12  h  over  the  100-km  width 
coastal  area  off  Peru  (Fig.  5e  and  f).  The  diurnal  amplitude 
along  the  Peruvian  coast  is  only  0.2-0. 5  ms-1,  along  with  a 
small  correlation  coefficient  about  0.2,  which  implies  signif¬ 
icant  uncertainties  in  the  model  diurnal  statistics.  The  merid¬ 
ional  variation  of  wind  speed  phase  in  this  area  has  been 
noticed  in  previous  studies.  For  example,  Munoz  (2008) 


used  model  simulations  to  show  a  decrease  in  meridional 
wind  speed  along  the  southern  coast  of  Peru  from  12:00  to 
18:00  FST  compared  to  the  corresponding  increase  over  ar¬ 
eas  off  central  coast  of  Chile.  Analysis  of  satellite  data  also 
suggest  that  the  timing  of  the  maximum  meridional  winds 
off  Peru  and  northern  Chile  is  generally  several  hours  earlier 
than  off  the  central  region  of  Chile  (e.g.,  Gille  et  al.,  2005). 
The  characteristics  of  the  diurnal  variation  of  surface  winds 
in  the  coastal  area  are  important,  as  they  reflect  the  interac¬ 
tion  among  the  lower  troposphere  circulation,  sea  breeze  and 
MBF  turbulence.  Given  the  high  uncertainties  regarding  the 
forecast- satellite  wind  phase  difference  and  the  meridional 
variation  of  the  wind  phase  shown  in  these  studies,  more  de¬ 
tailed  analyses  are  needed  to  produce  robust  diurnal  statistics 
of  both  model  and  observations  and  investigate  the  exact  na¬ 
ture  of  the  wind  diurnal  variation. 

This  regional  comparison  demonstrates  that  the  predicted 
mean  wind  and  cloud  fields  in  general  agree  with  the  satellite 
observations,  including  the  low-level  jet  along  the  west  coast 
of  Chile,  the  weak  wind  speed  area  present  near  the  subtrop¬ 
ical  high  center  and  in  the  coastal  area  near  20°  S,  the  overall 
FWP  distribution  pattern,  and  diurnal  phase  pattern  for  both 
clouds  and  wind  speeds.  The  predicted  diurnal  variance  of 
the  FWP  from  COAMPS  explains  10-25%  of  the  total  vari¬ 
ance,  a  percentage  smaller  than  that  (30-60%)  derived  from 
the  satellite  data.  Diurnal  variation  of  MBF  clouds  in  the 
southeast  Pacific  may  also  include  a  6-h  harmonic  which  is 
driven  by  the  so-called  upsidence  wave  as  discussed  by  Gar- 
reaud  and  Munoz  (2004)  and  O’Dell  et  al.  (2008).  This  vari¬ 
ability,  however,  is  not  included  in  the  current  analysis.  It 
is  discussed  in  detail  by  Rahn  and  Garreaud  (2010a)  using 
model  simulations  and  observations  in  VOCAFS-REx. 


4  20°  S  MBL  cross-section 

This  section  focuses  on  the  evaluation  of  the  COAMPS  MBL 
structure  on  the  15  km  grid  mesh  along  20°  S  using  the  RHB 
rawindsondes,  NCAR  C-130,  and  CIRPAS  Twin-Otter  mea¬ 
surements.  The  RHB  had  two  cruise  legs,  with  first  running 
from  6  October  to  3  November  and  the  second  from  1 1  to  29 
November.  No  data  were  collected  3-10  November  between 
the  two  periods  while  it  was  stationed  in  Arica,  Chile,  and 
only  the  data  collected  along  20°  S  are  used  here. 

Figure  6  compares  the  time-height  cross-section  of  po¬ 
tential  temperature  (i 9 )  vertical  profiles  of  the  lower  atmo¬ 
sphere  from  the  RHB  with  those  from  the  forecast.  This 
predicted  cross-section  is  derived  using  interpolation  based 
on  the  times  and  locations  of  the  rawindsondes.  The  com¬ 
parison  starts  at  23:24  UTC  25  October  2008  and  ends  at 
16:00  UTC  28  November  2008,  including  196  profiles  from 
each  dataset.  The  observed  nearly  well-mixed  MBL  structure 
capped  by  a  strong  inversion  is  well  predicted.  The  simulated 
jump  in  0  across  the  inversion  appears  to  be  comparable  to 
the  observed,  although  the  modeled  inversion  layer  is  thicker, 
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Fig.  6.  Time-height  cross-section  of  potential  temperature  ( 0 )  from 
the  COAMPS  forecasts  (left)  and  the  RHB  rawindsondes  (right). 
The  analysis  covers  the  period  of  25  October  through  29  November 
2008.  The  numbers  at  the  upper  axis  with  the  heavy  vertical  lines 
denote  the  longitudes  of  the  most  western  and  eastern  locations  of 
the  RHB  cruise  for  each  leg  during  this  period.  No  data  were  col¬ 
lected  3-10  November  between  the  two  cruise  legs. 


partly  due  to  the  limited  vertical  resolution  in  the  model.  For 
the  first  9  days,  the  variation  of  the  predicted  Zi  closely  fol¬ 
lowed  that  from  the  observations  except  that  it  lowered  to 
600  m  while  the  RHB  z,i  stays  above  1  km  near  the  coast.  At 
the  start  of  leg  2,  the  modeled  Zi  was  again  low  at  600  m, 
and  then  immediately  increased  to  1.4  km  and  varied  little 
afterwards.  In  contrast,  the  RHB  zi  stayed  at  ~  1.3  km  until 
1 8  November,  after  which  time  it  changed  considerably  and 
reached  a  maximum  of  2  km  at  20°  S  and  85°  W.  The  mod¬ 
eled  Zi  clearly  lacks  the  variability  comparable  to  the  RHB 
observations  in  leg  2.  The  modeled  Zi  is  generally  lower  than 
in  the  RHB  observations,  particularly  in  the  nearshore  area. 

During  VOCALS-REx,  two  different  dominating  synoptic 
conditions  are  identified  (Rahn  and  Garreaud,  2010b).  Dur¬ 
ing  October,  synoptic  disturbances  are  prominent,  strongly 
modulating  Zi  as  demonstrated  by  both  the  observations  and 
forecasts  as  shown  by  the  low  frequency  variability  in  the 
time-height  sections  of  potential  temperature  in  Fig.  6.  For 
this  period,  the  observed  variation  of  z,i  is  very  well  simu¬ 
lated.  In  contrast,  during  November,  synoptic  forcing  is  weak 
over  the  VOCALS-REx  region.  For  this  period,  COAMPS 
appears  to  have  difficulty  in  predicting  correct  Zi  near  20°  S 
and  85°  W.  A  possible  explanation  for  these  different  model 
behaviors  is  as  follows.  In  October,  tendency  of  Zi  may  be 
dominated  by  synoptic  disturbances  which  can  be  resolved 
by  the  model.  In  November  without  strong  synoptic  forcing, 
the  tendency  may  be  mostly  determined  by  the  turbulence 
or/and  shallow  convection  which  must  be  parameterized.  It 
is  possible  that  deficiencies  in  the  turbulence  parameteriza¬ 
tion  may  play  more  prominant  roles,  leading  to  more  errors 
in  Zi  during  November.  To  test  this  hypothesis,  one  needs  to 
analyze  the  budget  of  Zi  to  compare  the  relative  importance 
of  different  processes  in  the  Zi  budget. 


To  further  examine  the  MBL  vertical  structure,  we  com¬ 
posite  the  vertical  soundings  of  potential  temperature  0 , 
specific  humidity  qv ,  equivalent  potential  temperature  0e, 
and  the  zonal  and  meridional  wind  components  ( u  and 
v,  respectively)  according  to  their  distances  to  the  coast. 
Figure  7  compares  the  composite  soundings  centered  at 
3  locations:  “nearshore”  (71.64°  W,  19.22°  S),  “midway” 
(79.16°  W,  19.22°  S),  and  “offshore”  (84.5°  W,  20.4°  S). 
These  composites  are  produced  by  averaging  the  soundings 
within  one  degree  in  longitude  of  the  locations.  Most  of  the 
predicted  variables  closely  follow  the  observations  within  the 
MBL  and  above  the  inversions.  The  well-mixed  0  profiles  in 
the  lower  MBL  are  well  predicted;  the  inversion  strength  also 
compares  well  with  the  observations.  There  is  a  clear  trend 
of  increasing  moisture  above  the  inversion  toward  the  coast 
in  both  the  forecasts  and  the  RHB  observations.  A  weak¬ 
ness  in  the  predicted  MBL  structure  is  the  lower  Zi  at  both 
the  “nearshore”  and  “offshore”  locations;  particularly,  Zi  is 
underpredicted  nearshore  by  450  m.  It  should  be  mentioned 
that  the  41 -day  ensemble  averaged  Zi  at  85°  W  and  20°  S  is 
similar  to  the  C-130  soundings,  while  the  low  bias  of  the 
nearshore  Zi  remains  a  consistent  feature  of  the  forecasts, 
as  shown  in  next  section.  The  model  also  produced  higher 
moisture  just  above  the  inversion  and  below  2  km,  as  shown 
by  the  specific  humidity  ( qv )  comparison  in  Fig.  7. 

The  observed  equivalent  potential  temperature  ( 0e )  de¬ 
creases  with  height  below  the  inversion  top  (Fig.  7).  The  neg¬ 
ative  6  e  gradient  becomes  stronger  toward  the  open  ocean, 
suggesting  that  the  MBL  tends  to  be  more  conditionally  un¬ 
stable  offshore,  a  condition  favorable  for  shallow  cumulus 
development  (Stevens,  2007).  The  model  prediction  of  0e 
has  a  similar  feature  for  the  “offshore”  and  “nearshore”  lo¬ 
cations.  Due  to  the  higher  values  of  qv  just  above  the  inver¬ 
sion,  the  overall  change  in  the  modeled  0 e  across  the  MBL 
is  positive  at  the  “midway”  location,  in  contrast  to  the  RHB 
soundings. 

Both  the  model  and  the  RHB  observations  show  wind 
maxima  just  below  the  inversion  with  diminishing  values  up¬ 
ward  (Fig.  7).  The  low  wind  speed  above  the  MBL  is  asso¬ 
ciated  with  the  shallowness  typical  of  subtropical  high  pres¬ 
sure  systems.  The  meridional  wind  shear  across  the  inver¬ 
sion  is  consistent  with  the  zonal  temperature  gradient  within 
the  inversion,  introduced  by  the  east-west  upward  slope  of 
the  MBL  height  along  20°  S.  We  can  estimate  the  jump  in  v 
across  the  inversion  based  on  the  thermal  wind  relationship 
using  both  the  model  and  RHB  soundings.  First,  0  should 
remain  the  same  along  the  inversion  slope.  That  is, 

dzi(x)  =  0.  (1) 

Using  Eq.  (1)  and  the  model  soundings,  the  mean  zonal  tem¬ 
perature  gradient  within  the  inversion  can  be  approximated 
as, 


(d0)Zi  = 


dx  4- 
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Fig.  7.  Comparison  of  the  composite  COAMPS  and  RHB  soundings.  The  soundings  are  averaged  at  the  central  location  indicated  by  the 
pair  of  longitude  and  latitude  in  the  first  columns  in  the  figure.  The  digits  under  the  locations  are  the  total  number  of  soundings  used  in  each 
averaging.  The  variables  are  denoted  by  x-axis  labels.  Red  curves  denote  the  forecast  results;  and  the  black  the  RHB. 
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This  horizontal  temperature  gradient  gives  an  estimate  of 
the  mean  jump  in  v  using  the  thermal  wind  relationship: 


AUmodel  ~  —  •  (~ 

fO  \dx 


•1.2  x  10-5  -400 


A  Zi 


9.8 


-5  x  10“5  x  300 


-3  (ms  *), 


(3) 


which  is  similar  to  the  forecast  value  ~2-5ms_1  (Fig.  7). 
The  mean  jump  in  the  RHB  Aurhb  across  the  inversion  can 
also  be  estimated  in  a  similar  fashion.  That  is,  Aurhb  ~ 
— 4ms_1  for  a  12 K  increase  in  0  across  the  inversion  and 
a  700  m  increase  in  Zi  over  13  degrees  of  longitude.  This 
value  is  at  the  lower  end  of  the  observed  range  from  —4  to 
—7 ms-1  shown  in  Fig.  7,  and  is  consistent  with  the  esti¬ 
mate  by  Bretherton  et  al.  (2010)  using  multi-platform  obser¬ 
vations.  The  jumps  from  the  model,  however,  are  weaker 
than  those  from  the  observations  at  both  the  “offshore”  and 
“midway’  locations,  while  it  is  stronger  at  the  “nearshore”. 
This  is  consistent  with  an  underprediction  of  Zi  both  near  the 


coast  and  at  the  most  western  locations,  resulting  in  a  gentler 
slope  of  Zi  offshore  and  a  steeper  one  nearshore. 

CIRPAS  Twin-Otter  research  flights  provide  datasets  of 
coastal  MBF  structure  at  the  fixed  location  (20°  S,  72°  W). 
Figure  8  compares  mean  Zi  -normalized  boundary  layer  pro¬ 
files  collected  from  18  research  flights  with  the  model  fore¬ 
casts  at  the  corresponding  times.  We  first  note  that  the  mean 
predicted  MBF  height  is  751  m,  which  is  about  400  m  lower 
than  that  derived  from  the  Twin  Otter  as  shown  in  Fig.  8a, 
and  consistent  with  the  underprediction  of  zi  found  in  the 
nearshore  RHB  soundings  shown  in  Fig.  7.  The  mean  pre¬ 
dicted  temperature  and  moisture  within  the  MBF  in  gen¬ 
eral,  agrees  well  with  the  Twin  Otter  observations.  Both 
the  model  results  and  Twin  Otter  observations  show  appar¬ 
ent  wind  speed  and  directional  shear  across  the  inversion,  al¬ 
though  the  observed  changes  are  significantly  stronger  than 
the  forecast.  The  wind  direction  changes  from  southerly 
within  the  MBF  to  northerly  just  above  the  inversion  in  both 
the  observations  and  model  simulations.  This  wind  direc¬ 
tional  shear  is  also  implied  by  the  narrow  band  (about  500  km 
wide)  of  northerly  flow  above  the  MBF  as  shown  in  Fig.  2f 
near  20°  S,  72°  W.  Consequently,  this  wind  shear  is  associ¬ 
ated  with  the  meso-  or  large-scale  dynamics,  a  key  difference 
from  the  usual  wind  shear  across  the  inversion,  driven  by  tur¬ 
bulence  mixing.  This  wind  shear  may  introduce  significant 
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Avg.  COAMPS  vs.  Twin-Otter  (20S,  72W)  COAMPS  . 


Twin-Otter- 


0(K) 


Fig.  8.  Comparison  between  COAMPS  and  Twin  Otter  averaged  soundings  scaled  by  the  boundary  layer  heights,  (a)  0\  (b)  qv  and  10  xgc; 
(c)  wind  speed;  and  (d)  wind  direction. 


turbulence  mixing  within  the  entrainment  zone,  and  thus  has 
important  implications  on  the  MBL  turbulence  structure  and 
cloud-top  entrainment  rate  (Wang  et  al.,  2008).  It  should  be 
noted  that,  near  the  surface,  the  wind  direction  is  ~  175°  for 
the  Twin-Otter  and  ~210°  for  the  model.  Many  factors  may 
contribute  to  the  difference,  including  the  timing  and  strength 
of  the  sea  breeze,  movement  of  the  subtropical  high  pressure, 
impacts  of  steep  Andes  Mountains,  and  diurnal  variability  of 
the  stratocumulus.  An  examination  of  the  results  on  the  5 -km 
grid  mesh  indicates  that  the  41 -day  ensemble  mean  wind  di¬ 
rection  (not  shown  here)  is  203°,  a  slight  improvement  com¬ 
pared  with  that  on  the  15 -km  grid,  but  not  enough  to  correct 
the  overprediction  of  the  onshore  component. 

Because  of  the  models’  underprediction  of  Zi  compared  to 
observations,  the  simulated  0  just  above  the  MBL  is  lower 
than  that  of  the  Twin-Otter  observations.  The  air  is  also 
moister  than  the  observed  by  3gkg_1,  further  confirming 
the  early  result  that  the  air  just  above  the  inversion  is  too 
moist  near  the  coast  in  the  model  compared  with  the  RHB 
soundings  (Lig.  8b).  The  maximum  predicted  cloud  water 
(qc)  is  only  0.1  g  kg-1  compared  with  the  observed  value, 
0.3  g  kg-1.  The  low  values  of  qc  are  mainly  due  to  the  thin¬ 
ner  cloud  layer  caused  by  the  lower  Zi .  This  explains  the  low 
values  of  LWP  in  the  forecasts  near  the  coastal  area  shown  in 
Lig.  3a  and  Lig.  4. 

Longitude-height  cross-sections  of  the  simulated  41 -day 
time-averaged  MBL  structure  along  20°  S  are  shown  in 
Lig.  9.  Also  plotted  is  the  MBL  height  derived  from  all 
the  NCAR  C-130  20°  S  relative  humidity  and  radar  measure¬ 
ments  (C.  S.  Bretherton,  personal  communication,  2010).  In 
general,  the  predicted  Zi  increases  offshore  with  increased 
cloud  layer  depth.  It  agrees  with  the  observed  Zi  reasonably 
well  offshore  between  85°  W  to  75°  W,  particularly  between 
78°  W  and  75°  W.  The  predicted  zi,  however,  falls  sharply 
near  the  coast  to  the  east  of  75°  W,  consistent  with  the  com¬ 
parison  with  the  RHB  and  Twin-Otter  soundings  (Ligs.  7  and 
8.)  The  average  cloud  base  in  the  model  is  higher  than  the 
lifting  condensation  level  (LCL).  This  difference  (cloud  base 
level  minus  LCL)  sharply  increases  from  the  coast  to  about 
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Fig.  9.  COAMPS  MBL  variations  along  20°  S.  (a)  Longitude- 
height  cross-section  of  qc  (white  contours  in  gkg-1)  and  0  (color 
shading  in  K).  The  mean  predicted  MBL  height  is  denoted  by  the 
upper  branch  of  black  curve,  the  mean  LCL  by  the  lower  branch  of 
the  black  curve,  mean  cloud  base  by  the  black  dashed  curve.  The 
blue  line  is  the  MBL  height  fitted  by  C-130  data,  (b)  cross  sec¬ 
tions  of  qv  (white  contours  in  gkg-1)  and  0e  (color  shading  in  K); 
(c)  cross  sections  of  downward  LW  flux  (color  shadeing  in  Wm“2) 
and  upward  LW  flux  (white  contours);  and  (d)  predicted  longitudi¬ 
nal  variations  of  the  lower  troposphere  stability  ( A0  =  $3  km  —  $SST> 
solid  line  in  K),  the  time  averaged  vertical  motion  at  the  inversion 
height  (dashed  line  in  mm  s-1)  and  the  total  LW  divergence  across 
the  cloud  layer  (ALW,  dotted  in  Wm-2).  ALW  is  defined  as  the  dif¬ 
ference  between  the  total  LW  radiative  flux  just  above  clouds  and 
that  just  below  the  cloud  base. 


76°  W  and  then  increases  more  gently  toward  the  open  ocean, 
suggesting  that  the  MBL  becomes  decoupled  to  the  west  of 
76°  W.  That  is,  the  MBL  thermodynamic  profiles  depart  from 
being  well-mixed  and  become  more  stratified  due  to  stronger 
cloud-top  entrainment  (Bretherton  and  Wyant,  1997). 
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The  predicted  moisture  above  the  inversion  stays  constant 
at  about  2  gkg-1  offshore  (Fig.  9b);  it  increases  significantly 
to  the  east  of  76°  W.  Because  of  this  increase  in  moisture,  the 
downward  LW  flux  above  clouds  increases  from  230  Wm-2 
at  80°  W  to  300 Wm-2  at  72°  W,  as  shown  in  Fig.  9c.  The 
strong  gradient  in  the  downward  FW  flux  at  the  cloud  top 
indicates  large  LW  divergence,  producing  the  intense  LW 
radiative  cooling  that  is  a  main  driver  of  the  turbulence  in 
the  MBL.  The  minimum  upward  LW  flux  near  the  cloud  top 
reflects  the  effect  of  minimum  temperatures  at  these  levels. 
Owing  primarily  to  the  increase  in  the  moisture  above  the 
cloud  toward  the  coast,  the  LW  divergence  decreases,  re¬ 
sulting  in  the  reduced  total  LW  forcing  in  the  cloud  layer 
nearshore,  as  demonstrated  in  Fig.  9d.  This  trend  agrees 
with  the  RHB  and  NCAR  C-130  observations;  it  has  dy¬ 
namic  implications  on  the  MBL  turbulence  and  cloud  struc¬ 
ture.  Bretherton  et  al.  (2010)  presents  observational  evidence 
that  the  westward  decrease  in  moisture  above  clouds  con¬ 
tributes  to  the  westward  increase  in  the  LW  radiative  diver¬ 
gence  in  the  MBL,  and  consequently,  to  the  corresponding 
increase  in  the  turbulence  intensity  and  the  MBL  height. 

The  overall  increase  in  the  MBL  height  and  the  more  de¬ 
coupled  structure  offshore  are  both  consistent  with  an  off¬ 
shore  decrease  of  lower  troposphere  stability,  as  defined  here 
by  A0  =  Oikm  —  Osst  (Fig.  9d)  (Bretherton  and  Wyant,  1997; 
Klein,  1997).  The  41 -day  averaged  predicted  vertical  veloc¬ 
ity  at  the  inversion  height  level,  presented  in  Fig.  9d,  reaches 
a  minimum  (~  —  0.81cms_1)  at  72.3°  W  compared  to  the 
values  of  —0.3  to  —0.4  cm  s-1  to  the  west  of  75°  W.  This  sug¬ 
gests  that  the  unusually  strong  subsidence  in  the  model  fore¬ 
casts  is  likely  a  main  contributor  to  the  thinning  of  the  MBL 
along  the  coast.  Further  analysis  will  be  given  in  Sect.  7  on 
the  issue  of  the  predicted  low  MBL  heights  nearshore. 

In  Sect.  3,  the  cloud  diurnal  amplitude  field  was  shown  to 
be  spotty  compared  to  satellite  observations  (Fig.  4a).  Sim¬ 
ilar  spotty  cloud  fields  were  also  seen  in  the  real-time  fore¬ 
cast  presentations  during  VOCALS-REx  daily  weather  brief¬ 
ings.  There  is  strong  evidence  that  these  sporadic  cloud  fields 
are  likely  linked  to  the  so-called  grid-scale  convection  due 
to  the  lack  of  a  subgrid- scale  shallow  cumulus  convection 
parameterization.  These  sporadic  cloud  features  are  more 
clearly  seen  in  Fig.  10a,  which  displays  a  longitude-time 
cross-section  of  the  predicted  LWP  along  20°  S.  Apparent  are 
the  appearance  of  short-lived  and  thick  clouds  (LWP  ~500- 
700  g  m-2)  in  the  background  of  a  well-defined  diurnal  varia¬ 
tion  of  relatively  thin  clouds  (LWP  ~100-150gm-2).  Most 
of  the  spotty  and  thick  clouds  form  offshore,  where  0e  and 
qv  in  the  lower  MBL  increase  westward  and  decrease  with 
height  due  to  the  increase  in  SST  (Figs.  2a  and  9b).  The 
decoupled  MBL  offshore  also  tends  to  build  up  moisture  in 
the  subcloud  layer  by  inhibiting  upward  moisture  transport 
from  the  subcloud  layers.  These  are  favorable  conditions  for 
the  formation  of  shallow  cumulus  clouds.  Since  there  is  no 
shallow  cumulus  convection  parameterization  in  the  model, 
the  instability  must  then  be  released  by  grid-scale  convec- 


Fig.  10.  Examples  of  grid-scale  convection  in  the  COAMPS  sim¬ 
ulation.  (a)  Longitude-time  cross-section  of  the  COAMPS  LWP 
along  20°  S,  (b)  longitude-height  cross-section  of  the  instantaneous 
qc  (color  shading)  and  w (contours  starting  at  Ocms-1  with  an  in¬ 
terval  of  5cms_1)  along  20°  S  at  03:00  UTC  5  November  2008 
[denoted  by  the  white  bar  on  the  right  axis  in  a],  (c)  mean  profiles 
of  qv  and  qc  averaged  between  89°  W  and  80°  W,  and  (d)  mean 
profile  of  0. 

tion  in  the  MBL,  which  leads  to  the  formation  of  the  spotty 
and  relatively  thick  clouds.  A  snapshot  of  these  clouds  is 
provided  in  Fig.  10b,  where  grid-scale  updrafts  of  15  cm  s-1 
produce  clouds  with  1 .5  km  in  depth  and  70  km  in  width.  The 
overall  MBL  structure  averaged  between  90°  W  and  80°  W  at 
this  time  (Fig.  10c  and  d)  shows  a  common  shallow  cumu¬ 
lus  thermodynamic  structure  consisting  of  a  surface-based 
mixed  layer,  a  thin  stable  layer  at  the  cloud  base,  a  cloud 
layer  and  an  inversion  layer  at  the  cloud  top  (e.g.  Albrecht  et 
al.,  1995;  and  Stevens,  2007). 

The  grid- scale  convection  may  increase  the  diurnal  ampli¬ 
tude  if  it  occurs  in  the  early  morning  when  the  cloud  LWP 
is  likely  at  the  peak;  it  may  also  interrupt  the  natural  diurnal 
cycle  by  producing  additional  thick  clouds  in  the  afternoon. 
Both  scenarios  occur  as  the  episodic  convection  is  produced 
during  both  day  and  night  as  demonstrated  in  Fig.  10a.  For 
areas  where  the  convection  is  not  frequently  present  (e.g.,  to 
the  east  of  75°  W  and  between  19°  S  and  24°  S  from  Figs.  4a- 
d  and  10a),  the  forecast  diurnal  amplitude  and  correlation 
coefficient  appear  to  be  in  better  agreement  with  the  satellite 
data.  A  shallow  cumulus  parameterization  like  that  devel¬ 
oped  by  Park  and  Bretherton  (2009)  or  by  Golaz  et  al.  (2002) 
may  be  implemented  in  COAMPS  to  address  this  issue  with 
anomalous  grid-scale  shallow  convection. 

In  summary,  the  COAMPS  model  forecast  captures  the  ba¬ 
sic  large-scale  variability  of  the  MBL  along  20°  S,  includ¬ 
ing  the  general  increase  in  the  MBL  heights  toward  the  open 
ocean,  the  more  decoupled  MBL  structure  offshore  and  the 
directional  wind  shear  near  the  coast.  The  upward  slope  of 
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the  MBL  height  toward  the  west  is  consistent  with  the  de¬ 
crease  in  meridional  wind  speed  across  the  inversion  based 
on  the  thermal  wind  relationship.  The  analysis  of  the  instan¬ 
taneous  and  averaged  MBL  structure  suggest  that  the  spo¬ 
radic  features  in  the  cloud  field  (e.g.  Ligs.  4a  and  10a)  is 
likely  a  result  of  the  lack  of  shallow  cumulus  convection  pa¬ 
rameterization  in  the  model.  The  comparisons  with  RHB 
and  aircraft  soundings  consistently  show  a  relatively  large 
underprediction  of  the  MBL  heights  in  the  nearshore  area, 
which  may  adversely  impact  the  cloud  prediction.  We  show 
in  Sect.  7  that  this  underprediction  in  MBL  heights  near  the 
coast  by  the  model,  in  part,  result  from  the  coarse  horizontal 
and  vertical  resolution,  which  notably  affects  the  vertical  mo¬ 
tion  and  horizontal  advection  associated  with  the  steep  slope 
of  the  Andes. 


5  Surface  statistics 

This  section  focuses  on  the  evaluation  of  predicted  surface 
variables  with  the  observations  collected  during  the  RHB 
cruises  and  at  the  WHOI  buoy  moored  at  20°  S,  85°  W.  Lig- 
ure  11  displays  longitudinal  variations  of  surface  variables 
from  the  RHB  observations  and  model  forecast  data  along 
20°  S.  To  better  display  the  general  trend  of  the  variables, 
we  first  group  the  data  into  1 -degree  longitude  bins  and  then 
average  them  over  each  group.  The  temperature,  wind  speed 
and  water  vapor  mixing  ratio  are  measured  at  15-18  m  height 
on  the  RHB,  and  therefore  they  are  converted  to  the  10-m 
level  temperature  according  to  the  Monin-Obuhkov  Simi¬ 
larity  theory.  The  model  NCODA  SST  closely  follows  the 
RHB  SST  variation,  reaching  local  minimum  values  near 
78°  W  and  75°  W  and  increasing  toward  the  coast  as  well 
as  the  open  ocean.  It  is  also,  in  part,  responsible  for  the 
strong  lower  tropospheric  stability  between  78°  W  and  76°  W 
(Lig.  9d).  Reflecting  the  SST  variation,  both  the  predicted 
and  RHB  10-m  air  temperature  (7iom)  have  low  values  near 
78°  W  and  75°  W.  The  modeled  SST  and  7iom  have  biases 
of  0.22  and  —  0.67°  C,  respectively,  and  correlation  coeffi¬ 
cients  of  0.83  and  0.71,  respectively,  showing  good  agree¬ 
ment  with  the  observations.  The  predicted  10-m  specific  hu¬ 
midity  (qv io  m)  follows  the  RHB  near  the  coast  and  shows 
a  general  increasing  trend  with  the  increasing  SST  offshore 
(Lig.  11c).  However,  the  overall  correlation  is  low  (0.25)  and 
the  bias  is  0.44  g  kg-1.  The  predicted  wind  speed  increases 
from  4 ms-1  nearshore  to  7  ms-1  at  84°  W,  consistent  with 
the  RHB  observations.  Its  local  variation  is  also  comparable 
with  the  RHB  observations.  The  wind  speed  has  a  low  bias 
(—0.31  ms-1)  and  a  relatively  high  correlation  coefficient  of 
0.76. 

The  comparison  statistics  between  the  COAMPS  forecasts 
and  the  WHOI  buoy  are  presented  in  Table  1 .  The  tempera¬ 
ture  and  wind  statistics  are  in  general  very  similar  to  previous 
analyses  of  COAMPS  forecasts  over  the  area  along  the  US 
west  coast  (e.g.,  Doyle  et  al.,  2009).  The  shortwave  (SW) 


Fig.  11.  Comparisons  of  longitudinal  variations  at  20°  S  of  (a) 
SST,  (b)  T\ om>  (c)  <7i;10m>  and  (d)  10-m  wind  speed  between  the 
COAMPS  forecasts  and  the  RHB  observations.  All  the  plots  are 
produced  by  averaging  each  variable  into  1 -degree  longitude  bins 
using  the  corresponding  RHB  and  COAMPS  data.  The  overall 
statistics  (bias,  RMS  difference,  and  correlation  coefficient)  for  the 
COAMPS-RHB  comparisons  are  also  displayed  in  each  panel. 


downward  radiative  flux  from  the  forecasts  has  large  RMS 
error  compared  to  the  buoy  observations,  suggesting  that  the 
cloud-radiation  interaction  is  not  well  simulated.  The  long¬ 
wave  downwelling  flux  at  the  surface  correlates  poorly  with 
that  observed  at  the  buoy  (0.2  in  Table  1). 

To  better  understand  the  model  radiative  flux  errors,  we 
examine  scatter  plots  of  the  LW  fluxes  between  the  forecasts 
and  the  observations  collected  at  both  the  RHB  and  WHOI 
buoy  (Pig.  12).  Each  scatter  plot  displays  a  “square”  pattern 
with  very  similar  corresponding  values  at  the  four  corners. 
Presence  of  clouds  would  significantly  increase  the  down¬ 
ward  LW  flux.  Particularly,  when  the  cloud  layer  is  thick 
enough  to  radiate  as  a  blackbody,  the  thermal  radiation  by 
the  cloud  is  saturated  and  depends  only  on  the  cloud  layer 
temperature.  To  examine  the  effects  of  clouds,  we  separate 
the  “thick  clouds”  data  points  from  “clouds  free”  ones  using 
LWP  from  COAMPS  and  those  measured  by  the  RHB  ther¬ 
mal  infrared  radiometer.  The  “thick  cloud”  condition  (blue 
points)  is  defined  here  as  a  LWP  greater  than  30  gm“2  in 
both  the  observations  and  forecasts;  conversely,  the  “clouds 
free”  condition  (green  points)  is  defined  as  LWP  less  than 
1  gm“2  for  both  data  sets.  These  definitions  are  chosen  be¬ 
cause  they  successfully  stratify  the  data,  and  our  results  are 
not  particularly  sensitive  to  the  exact  cutoff  LWP  values.  As 
shown  in  Pig.  12a,  the  blue  cluster  at  the  upper-right  cor¬ 
ner  of  the  “square”  corresponds  to  the  condition  where  both 
the  observations  and  COAMPS  forecasts  are  designated  as 
“thick  clouds”,  while  the  green  cluster  at  the  lower-left  cor¬ 
ner  indicates  the  “clouds  free”  condition.  The  average  of  all 
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Table  1.  Statistical  comparisons  between  the  COAMPS  forecasts  and  the  WHOI  buoy  observations.  Correlation  coefficients  in  the  last  row 
do  not  have  any  units. 


T  (K) 

u  ( ms  1 ) 

v  (m  s  1 ) 

Speed 
(m  s->) 

LW  downward 
(W  m“2) 

SW  downward 
(W  m-2) 

Bias 

-0.43 

-0.51 

1.00 

0.35 

2.68 

-7.74 

RMS 

0.47 

1.54 

1.66 

1.35 

26.02 

188.68 

Corr. 

0.62 

0.52 

0.47 

0.62 

0.20 

0.84 

COAMPS-RHB  LW 


COAMPS-Buoy  LW 


Fig.  12.  Scatter  plots  of  surface  downward  LW  between  observa¬ 
tions  and  COAMPS  forecasts,  (a)  COAMPS  vs.  RHB;  and  (b) 
COAMPS  vs.  Buoy.  The  “clear”  condition  (green  points)  is  de¬ 
fined  here  as  LWP  less  than  1  g  m-2  from  both  the  CO  AMPS  and 
the  RHB;  The  “thick  cloud”  (blue  points)  is  defined  as  LWP  greater 
than  30  g  m-2  from  both  the  COAMPS  and  the  RHB. 


RHB  clear- sky  downward  LW  fluxes  of  315  W  m-2  is  close 
to  the  values  within  the  “clouds  free”  circle.  Since  the  major¬ 
ity  of  the  “blue”  and  “green”  data  points  are  located  along  the 
one-to-one  correlation  line,  the  LW  radiative  fluxes  for  both 
the  “clouds  free”  and  “thick  clouds”  conditions  are  therefore 
well  simulated  by  the  model.  These  clusters  of  points  ac¬ 
count  for  50%  of  all  data  points  (117  out  of  235  total).  Large 
radiative  flux  errors  occur  when  cloud  conditions  are  very 
different  between  the  forecast  and  observation  datasets.  For 
example,  the  points  at  the  upper-left  corner  with  values  cen¬ 
tered  at  (320,  380)  represent  the  “clouds  free”  condition  in 
the  observations  and  the  “thick  clouds”  in  the  COAMPS  fore¬ 
casts.  The  points  along  the  right-hand- side  vertical  line  of 
the  square  (i.e.,  RHB  =  380  Wm-2)  indicate  that  the  “thick 
clouds”  conditions  from  the  RHB  are  associated  with  thin 
cloud  or  partial  cloudiness  conditions  in  the  model.  This  LW 
comparison  raises  two  issues  related  to  the  cloud-radiation 
interaction  modeling.  One  is  the  realistic  prediction  of  qc  and 
partial  cloudiness;  the  other  is  the  physical  representation  of 
these  cloud  variables  in  the  radiation  calculation.  Both  issues 
need  to  be  further  addressed  to  improve  the  radiation  model¬ 
ing  in  the  future. 


6  Diurnal  variation  of  marine  boundary  layer  structure 

Section  2  presents  the  regional  cloud  variation  of  the  MBL 
diurnal  cycle  as  derived  from  both  the  model  forecasts  and 
satellite  data.  This  section  focuses  on  the  comparison  be¬ 
tween  the  forecast  and  observed  diurnal  cycle  of  MBL  struc¬ 
ture.  RHB  rawindsondes  were  launched  six  times  daily  at 
approximately  00:00,  04:00,  08:00,  12:00,  16:00,  and  20:00 
(UTC).  These  measurements  provide  adequate  sampling  of 
the  diurnal  variation  of  the  MBL  structure.  The  model  fore¬ 
cast  soundings  are  averaged  at  the  observation  times  to  pro¬ 
duce  the  mean  diurnal  cycles  for  both  the  forecasts  and  the 
RHB  observations.  To  reduce  the  effects  of  the  longitudinal 
variations  of  MBL  heights,  we  normalize  the  sounding  pro¬ 
files  using  the  MBL  heights  before  averaging. 

As  shown  in  Fig.  13,  both  the  MBL  heights  (zi)  calculated 
from  the  RHB  and  the  model  undergo  a  well-defined  diurnal 
cycle  with  the  maximum  in  the  morning  and  the  minimum 
in  the  afternoon.  This  diurnal  cycle  is  consistent  with  our 
general  understanding  that  the  daytime  solar  cloud  absorp¬ 
tion  tends  to  suppress  turbulent  intensity  and  reduce  cloud 
amount,  leading  to  weaker  cloud-top  entrainment  and  a  shal¬ 
lower  boundary  layer.  Both  the  forecast  and  the  observed  Zi 
reach  a  maximum  around  06:00  LST.  The  RHB  z,i  obtains  its 
minimum  at  15:00  LST,  4h  earlier  than  in  the  forecasts.  The 
RHB  Zi  is,  on  average,  140  m  higher  than  in  the  forecasts.  Its 
diurnal  amplitude  is  130  m  compared  to  only  60  m  for  that  of 
the  COAMPS  zi ,  which  is  consistent  with  our  earlier  result 
that  the  diurnal  amplitude  of  the  satellite  LWP  is  in  general 
larger  than  the  COAMPS  as  shown  in  Fig.  4. 

Figure  14  compares  the  MBL  structure  from  the  RHB  and 
the  forecasts.  To  better  present  the  results,  we  only  show 
the  profiles  at  06:54,  14:50  and  22:30  (LST).  The  general 
trend  of  the  MBL  warming  during  the  daytime  is  the  same 
in  both  the  forecasts  and  the  RHB  observations,  as  shown  by 
the  6  profiles  in  Fig.  14a  and  b.  The  MBL  temperature  in¬ 
creases  with  the  daytime  cloud-layer  solar  warming,  creating 
a  more  stratified  structure  at  14:50  LST.  In  the  evening,  the 
cloud  layer  cools  down  as  the  cloud  solar  warming  dimin¬ 
ishes  and  the  cloud-top  longwave  cooling  dominates  again. 
There  are  clear  differences  in  the  6  profiles  between  the 
observed  and  predicted  at  14:50  LST.  The  RHB  6  demon¬ 
strates  a  three-layered  structure  with  a  thin  stable  layer  at 
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Fig.  13.  Diurnal  variations  of  mean  MBL  heights  computed  from 
the  RHB  soundings  (black)  and  the  COAMPS  forecasts  (red). 


0.6  z/n,  separating  the  upper  stratocumulus  clouds  from 
the  surface-based  subcloud  layer.  The  COAMPS  0  profile 
has  a  smoothly  transitioned  two-layered  (cloud  and  subcloud 
layer)  structure.  These  differences  suggest  that  the  predicted 
daytime  MBL  has  a  less  decoupled  structure  than  the  obser¬ 
vations,  which  is  consistent  with  the  differences  in  the  cloud 
layer  warming  shown  between  the  forecast  soundings  (0.6  K 
between  07:00  and  14:50  LST)  and  the  RHB  (0.8  K)  shown 
in  Fig.  14. 

The  COAMPS  qv  is  larger  in  general  and  has  less  diurnal 
variability  than  the  RHB,  as  shown  in  Fig.  14c  and  d.  A  two- 
layer  structure  for  all  the  COAMPS  qv  soundings  is  clearly 
seen  with  the  subcloud  layer  being  nearly  well-mixed.  In 
contrast,  the  RHB  qv  profile  has  a  multi-layered  structure  at 
14:45  LST,  being  consistent  with  the  RHB  6  profile  at  the 
same  time.  The  maximum  diurnal  change  in  qv  from  the 
COAMPS  is  only  0.25  g  kg-1  compared  to  0.5  g  kg-1  from 
the  RHB  soundings. 

Associated  with  these  diurnal  changes  in  the  predicted  0 
and  qv ,  the  qc  and  TKE  also  undergo  similar  correspond¬ 
ing  changes,  as  shown  in  Fig.  14e  and  f.  The  qc  maximum 
is  0.35 gkg-1  at  06:54  LST,  and  decreases  to  0.15  gkg-1  at 
14:50  LST.  The  lower  cloud  layer  with  small  qc  represents 
shallow  cumulus  clouds  underneath  the  main  stratocumulus 
layer  between  0.6  and  1  z/n .  A  significant  reduction  in  TKE 
from  early  morning  to  afternoon  is  apparent,  particularly  in 
the  cloud  layer  (Fig.  14d).  The  minimum  values  of  the  TKE 
profiles  near  the  mid-MBL  indicate  that  the  cloud  layer  tends 
to  decouple  from  the  subcloud  layer. 


7  Sensitivity  of  MBL  structure  to  grid  resolution 

The  evaluation  of  the  COAMPS  forecasts  shows  that  the 
MBL  heights  close  to  the  coast  are  significantly  lower  than 
in  the  observations  (Figs.  7-9).  This  deficiency  appears  to  be 
quite  common  among  many  global  and  regional  models  as 


0  (K) 


0(K) 


Fig.  14.  Diurnal  variation  of  normalized  profiles  from  the  RHB  ob¬ 
servations  and  the  COAMPS  forecasts,  (a)  RHB  0\  (b)  COAMPS  0\ 
(c)  RHB  qv\  (d)  COAMPS  qv\  (e)  COAMPS  qc ;  and  (f)  COAMPS 
TKE.  The  numbers  in  (a)  denote  local  time  of  the  soundings  defined 
by  the  different  colors. 


discussed  by  Wyant  et  al.  (2010).  This  consistent  bias  among 
the  models  promotes  us  to  further  analyze  the  prediction  re¬ 
sults  at  different  horizontal  and  vertical  resolutions. 

We  first  compare  the  averaged  MBL  structure  on  three 
nested  grid  meshes  (45  km,  15  km,  and  5  km)  at  two  lo¬ 
cations:  nearshore  at  72°  W  and  20°  S,  and  offshore  at 
76°  W  and  20°  S.  By  increasing  the  horizontal  resolution 
from  45  km  to  5  km  ,  the  MBL  height  at  the  nearshore  lo¬ 
cation  increases  by  250  m  and  the  temperature  by  1  K  at 
72°  W  (Fig.  15).  Correspondingly,  qc  increases  by  50% 
to  0.15  g  kg-1  at  the  finest  grid  mesh  (5  km).  The  higher- 
resolution  grid  clearly  tends  to  improve  the  predicted  MBL 
structure  and  clouds  at  the  nearshore  location.  On  the  con¬ 
trary,  the  higher-resolution  grid  has  little  impact  on  the  MBL 
structure  at  the  offshore  location.  This  difference  implies 
that  the  improvement  for  the  nearshore  MBL  is  likely  caused 
by  changes  in  mesoscale  motions  due  to  the  high-resolution 
grid. 

Figure  15  further  shows  that  changes  in  downward  motion 
and  total  0  advection  between  the  coarse  and  fine  resolution 
grids  are  likely  responsible  for  these  structural  changes  at 
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Fig.  15.  Comparison  of  the  averaged  MBL  structure  and  mesoscale 
forcing  among  different  resolution  grids  and  locations.  Left  col¬ 
umn  of  panels  are  results  from  the  nearshore  location  at  72°  W  and 
20°  S;  the  right  column  the  offshore  location  at  16°  W  and  20°  S; 
the  horizontal  resolutions  (45  km,  15  km  and  5  km)  are  denoted  in 
the  top-right  panel,  the  letters  or  words  in  each  plot  denote  the  vari¬ 
ables. 


time  forecast.  The  48-h  sensitivity  simulation  is  initialized 
at  00:12  UTC  on  22  October  2008  and  its  time-averaged  re¬ 
sults  are  shown  in  Fig.  16.  Apparent  in  this  figure  is  that 
the  increase  in  vertical  resolution  leads  to  increases  of  MBL 
heights  for  all  three  grids  near  the  coast;  the  mean  MBL 
height  increases  the  most  for  the  5  km  grid,  with  differences 
of  ~300m.  The  increased  vertical  resolution  has  little  im¬ 
pact  on  the  MBL  structure  offshore  (not  shown  here).  As¬ 
sociated  with  the  increase  in  the  MBL  height,  the  maximum 
qc  increases  for  the  15  km  and  5  km  grid  due  to  the  cloud 
thickening.  Figure  16  also  makes  an  important  point  that 
the  high-vertical  resolution  produces  weaker  subsidence  in 
the  MBL  than  the  control  forecast  for  all  three  grids.  For 
the  finest  (5  km)  grid  mesh,  the  high-resolution  simulation 
gives  the  weakest  subsidence  and  produces  the  highest  MBL 
height.  There  is  therefore  a  clear  correlation  between  the  in¬ 
creased  MBL  heights  and  the  weakened  subsidence  within 
and  just  above  the  MBL  produced  by  the  increased  vertical 
and  horizontal  resolution  meshes. 

The  above  results  suggest  that  the  increases  in  both  the 
horizontal  and  vertical  resolutions  lead  to  the  improvements 
in  the  simulated  MBL  heights  and  clouds  near  the  coast,  and 
these  improvements  may  come  from  better  defined  coastal 
circulations.  Increasing  the  horizontal  resolution  better  re¬ 
solves  the  topography  and  coastal  flows.  The  increase  in 
the  vertical  resolution  likely  reduces  numerical  errors  asso¬ 
ciated  with  gradient  calculations  near  a  steep  terrain.  This 
is  consistent  with  some  previous  studies  demonstrating  that 
numerical  treatment  of  mountain  slopes  is  critical  for  a  real¬ 
istic  simulation  of  local  mesoscale  circulations  (e.g.  Zangl, 
2002).  However,  the  details  about  the  dependence  of  coastal 
mesoscale  circulations  on  the  model  resolutions  and  numeri¬ 
cal  techniques  are  beyond  the  scope  of  this  study. 


72°  W  and  20°  S.  The  finest-resolution  mesh  (5  km)  results 
in  the  weakest  subsidence  at  the  MBL  height,  leading  to  the 
highest  MBL  height  among  the  three  grid  meshes.  The  to¬ 
tal  3 -dimensional  0  advection  in  the  45  km  grid  mesh  gives 
the  largest  cooling  rate  near  the  surface;  that  is  double  that 
from  the  other  grid  meshes.  The  total  advection  near  the 
surface  mainly  represents  just  the  horizontal  component,  be¬ 
cause  the  vertical  advection  is  very  small.  Consequently,  this 
large  horizontal  advective  cooling  and  the  strong  wind  di¬ 
vergence  tend  to  strengthen  the  cloud-top  inversion  and  limit 
the  MBL  growth,  resulting  in  the  lowest  MBL  height  in  the 
45  km  grid  mesh.  It  should  be  noticed  that  the  subsidence  and 
the  total  advection  heating  are  very  similar  among  the  three 
grid  meshes  at  76°  W  and  20°  S,  which  is  consistent  with  the 
similar  MBL  heights  predicted  at  the  offshore  location. 

Effects  of  vertical  resolution  on  the  forecast  MBL  are  also 
evaluated  by  conducting  another  model  simulation  with  an 
enhanced  vertical  resolution  of  95  levels  (50  levels  below 
2  km)  compared  to  the  45  (23  levels  below  2  km)  in  the  real- 


8  Summary 

COAMPS  was  employed  to  provide  twice-daily  48  h  fore¬ 
casts  during  VOCALS-REx  in  October  and  November  2008. 
The  model  was  configured  with  three  nested  grid  meshes 
(45  km,  15  km  and  5  km)  and  45  vertical  levels  to  cover  the 
large  experiment  area  with  adequate  resolution.  The  MBL 
structure  from  the  real-time  forecasts  are  evaluated  and  vali¬ 
dated,  with  an  emphasis  on  results  from  the  15  km  grid  mesh, 
against  observations  taken  on  a  number  of  platforms,  includ¬ 
ing  aircraft,  ship,  buoy  and  satellites. 

The  regional  distributions  of  the  41 -day  averages  of  pre¬ 
dicted  LWP  and  wind  speed  are  generally  in  agreement  with 
those  observed  by  satellite,  although  the  modeled  LWP  is 
slightly  smaller  and  wind  speed  larger  in  the  VOCALS-REx 
region.  The  overall  location,  direction  and  strength  of  the 
southerly  coastal  jet  along  the  Chilean  coast  are  correctly 
predicted.  The  harmonic  regression  analysis  shows  that  the 
phase  of  diurnal  variation  of  the  MBL  clouds  and  the  surface 
wind  speed  are  in  general  consistent  with  the  satellite  data 
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Fig.  16.  Comparison  of  the  48  h-averaged  MBL  profiles  between  the  control  (45  levels)  and  the  high- vertical  resolution  (95  levels)  simulation 
at  72°  W  and  20°  S.  Upper  panels:  0\  and  middle  panels:  qc ;  and  the  lower:  w.  The  control  forecast  is  denoted  by  black  curves;  and  the  high 
resolution  by  the  red.  The  horizontal  resolution  is  45  km  for  the  left  panels,  15  km  for  the  central,  and  5  km  for  the  right. 


except  for  the  coastal  area  off  Peru.  The  diurnal  amplitude  of 
the  predicted  LWP  exhibits  a  spotty  spatial  distribution  com¬ 
pared  to  the  smooth  distribution  of  the  satellite  data.  We  have 
shown  that  this  spotty  distribution  is  caused  by  grid-shallow 
convection  likely  due  to  a  lack  of  a  shallow  cumulus  convec¬ 
tion  parameterization  in  the  model. 

The  forecast-observation  comparisons  along  20°  S  demon¬ 
strate  overall  consistency  in  some  basic  MBL  characteris¬ 
tics,  including  a  general  upward  slope  in  the  MBL  height, 
a  more  decoupled  boundary  layer  structure  toward  the  west, 
and  maximum  wind  speeds  near  the  top  of  the  MBL.  Both 
the  forecasts  and  the  observations  show  a  trend  of  decreas¬ 
ing  moisture  above  the  MBL  away  from  the  coast,  result¬ 
ing  in  stronger  radiative  cooling  near  the  cloud  tops  offshore. 
The  more  decoupled  MBL  structure  with  upward  slope  in  the 
MBL  height  toward  west  is  primarily  driven  by  the  westward 
decrease  in  the  lower  troposphere  stability,  which  enhances 
the  cloud-top  entrainment.  This  is  also  consistent  with  the  in¬ 
creased  LW  radiative  forcing  offshore.  The  model  analyzed 
SST,  predicted  10-m  air  temperature  and  surface  wind  speed 
yield  realistic  zonal  variations  along  20°  S  compared  with 
ship-based  RHB  observations.  The  model  consistently  pre¬ 
dicted  a  shallower  MBL  (~700m)  than  what  was  observed 
(~1 100  m)  in  the  coastal  region.  This  bias  leads  to  an  overall 
low  cloud  water  mixing  ratio  in  the  model  forecasts  com¬ 
pared  to  the  observations.  The  overall  statistical  errors  of 
the  surface  wind  speed  and  temperature  are  consistent  with 


previous  analyses,  with  the  wind  speed  bias  (<0.35 ms  l) 
and  RMS  (1.3  ms-1)  and  the  temperature  bias  (— 0.7  K)  and 
RMS  (0.5  K). 

The  simulated  LW  on  both  the  “thick  cloud”  and  “clouds 
free”  conditions  produce  good  agreements  with  the  observed 
values,  but  it  suffers  when  the  prediction  of  clouds  is  in  error, 
suggested  by  the  “square”  patterns  on  the  model-observation 
LW  scatter  plots.  This  clearly  shows  the  importance  of  the 
accurate  prediction  of  cloud  water  mass  and  fractional  cov¬ 
erage  in  the  determination  of  LW  radiative  fluxes. 

COAMPS  also  predicted  a  well-defined  diurnal  variation 
in  the  temporally  and  spatially  averaged  MBL  structure, 
which  is  generally  consistent  with  the  RHB  observed  n ,  6 
and  qv .  The  predicted  mean  structure,  however,  has  smaller 
diurnal  changes  in  these  variables  compared  to  the  observa¬ 
tions;  this  feature  is  also  reflected  in  the  smaller  diurnal  am¬ 
plitude  in  the  regional  LWP  diurnal  harmonic  analysis. 

To  further  understand  the  underprediction  of  MBL  heights 
near  the  coast,  we  examined  the  impacts  of  the  model  reso¬ 
lution.  It  was  found  that  higher  resolutions  in  both  the  hor¬ 
izontal  and  vertical  directions  indeed  have  noticeable  pos¬ 
itive  impacts  on  the  nearshore  MBL  heights,  while  having 
little  influence  offshore.  The  simulations  with  the  5 -km  grid 
mesh  (45  vertical  levels)  predicted  higher  and  more  realis¬ 
tic  MBL  height  (~900m)  than  those  with  the  15  km  mesh 
(~700  m)  or  the  45  km  mesh  (~600  m).  Enhancing  the  verti¬ 
cal  resolution  to  95  levels  further  increased  the  MBL  height 
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to  1100m  for  the  5  km  grid.  These  improvements  are  likely 
due  to  the  changes  to  the  mesoscale  flows  such  as  the  reduced 
coastal  downward  motion  and  weakened  cold  advection  near 
the  surface  produced  by  the  high-resolution  forecasts.  Given 
the  generality  of  this  bias  among  regional  and  global  models, 
and  the  implications  of  the  resolution  study  in  this  work,  ef¬ 
forts  should  be  directed  to  improve  simulations  of  mesoscale 
flows  under  the  influence  of  steep  terrains. 
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